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I .  INTRODUCTION 


One  of  the  vexing  problems  in  the  implementation  of  multidimen¬ 
sional  computer  codes  for  interior  ballistics  applications  is  the 
proper  handling  of  the  flow  singularity  at  the  conjunction  of  the  gun 
tube  wall  and  the  base  of  a  moving  projectile.  The  value  of  an  incor¬ 
rectly  determined  flow  variable  at  this  point  is  likely  to  be  propagated 
into  tie  interior  of  the  flow,  possibly  leading  to  disastrous  conse¬ 
quences  for  the  rest  of  the  flow  field.  In  the  past,  the  comer  point 
singularity  was  treated  by  ignoring  it,  i.e.,  using  a  coarse  mesh  with 
a  computational  stencil  which  does  not  include  the  comer  point,  or  by 
assuring  that  the  comer  is  multivalued,  i.e.,  it  is  both  at  rest  and 
moving  with  the  projectile  velocity.  Neither  of  these  approaches  is 
completely  satisfactory.  The  nature  of  the  flow  immediately  adjacent 
to  the  comer  should  suggest  the  correct  boundary  condition  at  the 
comer  point  itself.  Consequently,  we  apply  a  computational  stencil 
which  excludes  the  comer  point  and  examine  the  nature  of  and  changes 
in  the  flow  field  as  we  approach  the  sliding  comer  point  by  a  series 
of  grid  refinements. 


II.  THE  MODEL  PROBLEM 

Consider  a  smooth  bore  gun  tube  at  an  early  time  in  the  ballistic 
cycle.  The  projectile,  propelled  by  the  high  pressure  combustion  gases, 
proceeds  down  the  gun  tube.  For  the  case  at  hand,  the  obturation  is 
perfect,  blowby  is  not  allowed.  The  flow  field  behind  the  accelerating 
projectile  consists  of  a  core  flow,  boundary  layers  on  the  gun  tube 
walls  and  projectile  base,  and  the  comer  flow. 

Looking  at  a  slab  cut  out  of  this  tube  along  any  diameter,  we  have 
abstracted  the  following  configuration:  a  rectangular  region  bounded  on 
the  right  by  a  moving  projectile,  below  by  a  smooth  adiabatic  impermeable 
wall,  above  by  the  axis  of  symmetry  and  at  the  left  by  the  inflowing  gas 
with  a  given  velocity  distribution. 

The  flow  is  assumed  to  be  incompressible,  viscid  and  steady  through¬ 
out  the  region  of  interest.  From  the  Navier-Stokes  equations,  and  the 
definition  of  vorticity 
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one  obtains  the  vorticity  transport  equation 
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With  the  definition  of  the  stream  function 


ly  "  u*  lx  *  ”V  (3) 

equation  (1)  transforms  into  the  Poisson  equation 

V2*  -  C  •  (4) 

Introducing  a  characteristic  length,  L,  the  tube  radius,  and  a  time 
defined  as  the  advective  time  constant  (L/uo),  with  uQ  a  characteristic 

constant  projectile  speed,  the  equation  may  be  brought  into  the  following 
nondimensional  form 

;t  -  -Re  V  *(v  O*  ?2C  (5) 
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where  e  is  the  vorticity  nondimens ionali zed  by  (uq/L),  Re  is  the  Reynolds 

number,  v  the  vector  velocity.  Subscripts  denote  partial  differentiation. 
Since  we  are  looking  at  the  steady  case  only,  ■  0.  A  moving  coordi¬ 
nate  system  is  now  introduced  which  brings  the  projectile  to  rest.  One 
of  the  advantages  which  accrues  is  that  one  does  not  have  to  deal  with 
partial  cells  in  the  course  of  the  calculation. 


III.  THE  BOUNDARY  CONDITIONS 

Along  the  inflow  boundary,  Poiseuille  flow  and  thus  the  values  of 
the  stream  and  vorticity  functions  are  specified.  The  upper  boundary  is 
the  axis  of  the  tube,  a  line  of  symmetry  at  which  c  =  0  and  ♦  is  deter¬ 
mined  from  u  *  -P-  . 

max  3y 

In  the  projectile  fixed  coordinate  system  the  lower  boundary  moves  to 
the  left  at  a  velocity  of  u  ■  -1  and  ♦  ,  the  stream  function,  is  set  equal 
to  zero.  The  wall  vorticity  is  calculated  from  Woods'  equation; 

see  reference  1  for  details.  The  basic  idea  is  to  expand  the  stream 


T.  Roache,  P.J.,  Computational  Fluid  Dynamics,  Revised  Printing, 
Hermosa  Publishers,  Albuquerque,  New  Mexico,  1976. 
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function  in  a  Taylor  series  around  a  point  one  aesh  width  away  from  the 
wall  retaining  terms  up  to  the  third  order.  Using  the  definition  of 
vorticity  and  applying  the  no  slip  condition  on  the  surface  one  arrives 
at 


Ay2  " 

where  the  subscripts  w+1  and  w  denote  a  point  one  mesh  width  away  from 
the  wall  and  the  wall  point,  respectively.  A  similar  expression  holds 
along  the  piston  face. 

The  boundary  conditions  on  the  projectile  were  as  follows:  points 
along  the  projectile  base  have  zero  velocity  components;  thus  the  stream 
function  is  zero. 


IV.  TOE  ALGORITHM 

A  semidirect  solution  technique  was  used,  motivated  by  the  effi¬ 
ciency  of  these  methods.  As  in  (2),  the  equations  are  first  linearized, 
then  solved  directly;  subsequently  the  linearization  is  updated  and  the 
iteration  is  repeated.  Unlike  the  '"‘'f.hod  in  (1),  the  variable  coeffi¬ 
cient  advection  terms  are  included  in  the  direct  linear  solution  by  the 
use  of  marching  methods  (3) . 

The  marching  method  is  unstable  for  elliptic  partial  differential 
equations  as  the  mesh  is  refined,  i.e.,  as  AA->0  where  At  is  the  mesh 
size.  For  finite  mesh  widths,  however,  usable  results  may  be  obtained. 

Converged  answers  are  typically  obtained  in  less  than  10  iterations. 
The  total  computational  time  on  a  CDC  6600  is  of  the  order  of  5-10 
msec/cell.  No  empirical  or  semi-empirical  factors  (such  as  time  steps, 
under-relaxation  factors,  etc.)  need  to  be  determined. 


V.  TOE  CALCULATIONAL  SETUP 

In  a  planar  geometry,  a  series  of  calculations  were  carried  out 
with  the  number  of  grid  points  varying  from  11  X  H  to  41  X  41.  While 
the  radial  dimension  was  one  in  all  cases,  the  axial  dimension  was  in¬ 
creased  from  1,  as  the  number  of  cells  were  increased,  such  that  the 
cell  aspect  ratio,  axial  to  radial  dimension  varied  from  1:1  to  4:1. 

T.  Roache,  P.J.,  and  Ellis,  M.A. ,  "The  BID  Method  for  the  Steady-State 
Navier-Stokes  Equations",  Computers  and  Fluids,  3,  305-320  (1975). 

3.  Roache,  P.J.,  "Marching  Methods  for  Elliptic  Problems:  Part  1," 

Numerical  Heat  Transfer  1,  1-25  (1978).  "Part  2",  1,  163-181  (1978). 
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It  was  unfortunate  that  stability  requirements  of  the  algorithm  necessi¬ 
tated  the  use  of  these  cell  aspect  ratios  which  made  direct  comparison 
among  the  several  runs  difficult.  In  all  cases  the  wall  was  accelerated 
to  the  left  at  a  velocity  of  -1.  A  very  gentle  pressure  gradient,  of  the 
order  of  It,  was  imposed  in  the  axial  direction. 

To  start  the  calculation,  an  upstream  Poiseuille  flow  profile,  of 
2  3 

the  form  *  y-l.Sy  ♦  0.5y  ,  where  ^  is  the  stream  function  and  y  the 
radial  distance,  was  prescribed.  Analogous  prescriptions  were  made  for 
the  vorticity.  To  initialise,  the  Reynolds  number  was  set  to  zero,  then 
in  increments  of  10  it  was  raised  to  SO  in  the  course  of  the  calcula¬ 
tion  for  a  given  computational  mesh. 

Several  checks  monitored  the  progress  of  the  calculation.  If  the 
calculation  failed  to  converge  in  21  iterations,  the  calculation  was 
broken  off,  and  an  error  message  printed.  In  addition,  the  results  of 
two  convergence  checks  were  printed  out. 


VI.  DESCRIPTION  OF  THE  RESULTS 

Based  on  experimental  observations  and  classical  analysis,  one 
expects  the  Poiseuille  flow  to  persist  in  the  tube  coupled  with  a  region 
of  intense  vortex  generation  at  the  comer  itself.  This  expectation 
was  borne  out  by  the  numerical  calculations  which  we  now  proceed  to 
describe. 

The  results  are  best  illustrated  in  terms  of  plots  of  the  flow 
variables.  The  first  figure  shows  the  velocity  profile  at  the  position 
0.6  from  the  upstream  (inflow)  boundary  for  Re  »  30.  It  is  a  Poiseuille 
profile  in  the  transformed  coordinates. 

Next  the  contour  plot  of  the  u  velocity  component  gives  an  idea  of 
the  gradients  in  the  flow  velocity  in  the  computational  regime.  The 
contour  lines  are  spaced  at  0.2.  An  even  better  overall  view  can  be 
obtained  from  the  vector  velocity  plot.  Figure  3.  Note  that  the  flow 
(in  the  projectile  fixed  coordinates)  is  parallel  to  the  tube  axis, 
then  along  the  face  of  the  projectile,  finally  turning  around  and  exiting. 
Since  it  is  incompressible,  the  flow  in  and  out  of  the  region  of  interest 
must  balance  exactly. 

A  significant  insight  into  the  nature  of  the  flow  singularity  in 
the  comer,  of  course,  is  furnished  by  the  vorticity  plots.  Figures  4-7. 
The  contour  lines  are  spaced  at  2.0.  While  the  contour  plot  gives  a 
good  overall  idea  of  the  vorticity  distribution,  more  detail  is  furnished 
by  vorticity  plotted  along  the  wall  as  well  as  vorticity  along  the 
projectile  face. 
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In  both  instances  we  see  that  near  the  corner  singularity  the 
vorticity  rises  sharply.  This  is  readily  seen  in  the  three-dimensional 
plot,  Figure  7.  Also,  as  the  calculations  showed,  the  level  of  vortic¬ 
ity  is  increased  Minimally  as  the  Reynolds  number  was  raised  from  10  to 
50,  but  by  up  to  20%  in  going  to  a  31  X  31  mesh  from  an  11  X  11  (see 
Figures  6  and  8).  Away  from  the  comer  in  the  axial  direction,  practically 
there  was  no  variation  in  the  level  as  is  evident  in  Figure  6,  for 
example. 

The  results  further  show  that  the  lower  wall  is  contributing  the 
major  share  of  the  vorticity  generation.  Indeed,  if  one  considers  the 
ratio  of  the  vorticity  generated  by  the  lower  wall  to  that  generated 
along  the  projectile  base,  say  one  mesh  width  away  from  the  comer,  for 
a  mesh  of  31  X  31  at  Re  ■  10  the  ratio  is  1.8.  Similar  ratios  are  ob¬ 
tained  for  the  other  calculated  meshes. 

The  coarseness  of  the  mesh,  as  well  as  the  cell  aspect  ratio,  had  a 
definite  effect  on  the  maximum  level  of  vorticity.  As  the  number  of 
cells  was  increased,  the  cells  were  also  stretched;  the  resulting  mesh 
refinement  moved  the  cell  centers  further  away  from  the  comer  in  the 
axial  direction. 

One  should  note,  however,  that  these  results  must  be  interpreted 
carefully,  since  in  these  preliminary  runs  mesh  independence  of  the 
vorticity  was  not  obtained. 


VII.  CONCLUSIONS  AND  OUTLOOK 

At  this  point  our  conclusions  are  preliminary  in  nature.  We  found, 
through  mesh  refinement  near  the  confluence  of  the  projectile  and  a 
moving  wall, that  the  calculated  flow  there  conforms  to  qualitative 
estimates. 

As  one  proceeds  radially  from  the  comer,  the  velocity,  in  the 
transformed  coordinates,  is  seen  to  decrease.  This  is  true  in  the 
neighborhood  of  the  sector  formed  by  the  projectile  and  the  cylinder 
wall . 


The  vorticity  function  exhibits  a  sharp  rise  as  one  approaches  the 
comer  along  the  wall  in  the  axial  direction,  and  a  sharp  drop  when  the 
approach  is  from  the  radial  direction,  reminding  one  of  a  sink  type  flow 
behavior. 
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VELOCITY 


Figure  1.  Velocity  Profile  at  a  Position  0.6  from  the  Inflow 
Boundary  at  Re>30. 
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Figure  7.-  Vorticity  Surface  Plot  in  the  Computational  Domain. 
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